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The temperature (T)-shift protcol of aging in the 3 dimensional (3D) Edwards- Anderson (EA) 
spin-glass (SG) model is studied through the out-of-phase component of the ac susceptibility 
simulated by the Monte Carlo method. For processes with a small magnitude of the T-shift, 
AT, the memory imprinted before the T-shift is preserved under the T-change and the SG 
short-range order continuously grows after the T-shift, which we call the cumulative memory 
scenario. For a negative T-shift process with a large AT the deviation from the cumulative 
memory scenario has been observed for the first time in the numerical simulation. We attribute 
the phenomenon to the 'chaos effect' which, we argue, is qualitatively different from the so-called 
rejuvenation effect observed just after the T-shift. 

KEYWORDS: spin glass, slow dynamics, aging, droplet theory 
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§1. Introduction 

Recently- in studies on slow dynamics in spin 
glasses pel© the apparently contradictory phenomena, 
i.e., rejuvenation (or chaos) and memory effects in _a 
ing dynamics, have been intensively investigated!! 
In fact the phenomena were already observed in the 
early stage of study on aging in spin, glasses by the so- 
called temperature-cycling protocoLtf In the protocol 
we quench a spin glass to a temperature, say T\, be- 
low the spin-glass (SG) transition temperature T c from 
above it and let the system equilibrate (or age) for a pe- 
riod of t w i. Subsequently, we change the temperature 
to T 2 (< T\), age the system for a period of t w2 , and 
then the temperature is turned back to T\ . For a certain 
range of the parameters Ti, T%, i w i and i w2 some quanti- 
ties such as the out-of-phase component of ac suscepti- 
bility, y/'(o;;t), exhibit the following behavior. Just after 
the first T-shift y"(cL>;t) behaves as if the system were 
quenched to T 2 directly from above T c , or it looks hav- 
ing forgotten the aging at T\ before the T-shift. This is 
called the rejuvenation (or chaos) effect. However, after 
the temperature is turned back to Ti, x"(^;i) observed 
is the one we expect as a simple extension of x"(io;t) 
already aged by t w i at T\. Thus the system definitely 
preserves the memory of the previous aging at T\ , while 
it apparently exhibits the rejuvenation behavior, in the 
aging process at T 2 . The proper understanding of such a 
peculiar phenomenon is believed to shed light not only on 
the mechanism behind the aging dynamics but also the 
nature of the SG phase itself. Furthermore it will pro- 
vide us with powerful concepts to understand the glassy 
dynamics in various related systems such as orientational 
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glasses^ polymers^ and interacting nanoparticles sys- 
tems.B' 

By_tJjia pal-space interpretation, or by the droplet pic- 
t ur e,E2r E a c wh j ch we have been adopting in our recent 
studies,E3 , EJ't£lE3' the SG order is considered to grow up 
slowly in aging processes. In particular, we have demon- 
strated that the SG coherence length, which we regards 
as the mean size of SG domains developed in aging, con- 
tinuously grows even under the T-shift processE3' This 
we call the cumulative memory effect. We have further 
extended this characteristics to a scenario that the SG 
short-range order ever grows continuously with growth 
rates sensitively dependent on the_temperature so long 
as the system is in the SG phaselLa) Let us call this the 
cumulative memory scenario, and denote the mean size 
of SG domains as R T ^{t), where T[t] symbolically rep- 
resents the temperature changes that the system has ex- 
perienced up to time t from the first quench to the SG 
phase at t — 0. It has been demonstrated that the time 
evolution of zero-field-cooled (ZFC) magnetizations ob- 
served in various schedules of temperatuxe changes but 
with a common initial quench conditionliZP are well de- 
scribed by a unique function of Rx\t](t).cB 

The purpose of the present work is to numerically ex- 
plore to what extent the simulated data of x"(o;;i) in 
the T-shift protocol of the 3D Ising EA SG model are 
compatible or incompatible with the cumulative memory 
scenario, and with the rejuvenation (chaos) and memory 
effects mentioned above. For this purpose we certainly 
need a few more length and time scales than Rx[t}(t). 
One is Rx\t](t) at t = t w i, i.e., the mean domain size 
grown in the isothermal aging after quench to the SG 
phase, which is denoted as i?Ti(^wi)- After the temper- 
ature is shifted to T 2 spin configurations within each do- 
mains, which were in local equilibrium of Ti just before 
the T-shift, gradually become in local equilibrium of T 2 . 
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By the word 'gradually' we mean that the change asso- 
ciates with slowly growing droplets (or subdomains) of a 
mean size Lt 2 {t) (< i?Ti(iwi)) with r = t — t w i- Here 
we call this the droplets-in- domain scenario (previously 
called the quasi- domains-within- domains pictured) . 

A key quantity of the present work isJJne time scale re- 
quired for L T2 (t) to catch up Rt x (i w i ) » At time scales 
larger than this one, denoted as behavior of physi- 
cal quantities such as x"(Lj;i) cannot be distinguished, 
within the accuracy of measurement, from the corre- 
sponding behavior in the isothermal aging at T2. In other 
words, the system merges to a T2-isothermal aging state 
at the time scale of t^ after the T-shift. The latter is 
regarded as the effective waiting time of the i w i -aging at 
2~i reread as an T2-isothermal aging. If the cumulative 
memory scenario holds, is the time required for the 
SG coherence to grow in the T2-isothermal aging up to 
-RTi(iwi), i-e., 



Lt 2 (Ci ) = Rt 2 (f ■ 
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(1.1) 



is expected to hold. Here the first equality simply indi- 
cates that Tt(^) has the same functional form as that of 
Rrit) since the both growth processes are governed by 
common thermally-activated dynamics. 

In the present work we have extensively examined 
T-shift processes in the 3D Gaussian EA model with 
T c ~ O.OSJlMf through the ac susceptibility x"(w;t) sim- 
ulated by the standard Monte Carlo (MC) simulation. 
Here J is the variance of the interactions. The tempera- 
ture range investigated is T ~ [0.4, 0.7] in unit of J, and 
the time range is up to 10 5 MC steps. One of the results 
we have found is that in negative (positive) T-shift pro- 
tocols with Ti — » («— ) T2 Eq.(l.l) (the same equation 
but with the suffix 1, 2 interchanged) holds well when 
AT = T± — T% = 0.1. This confirms the cumulative na- 
ture of aging in both negative and positive T-shift proto- 
cols with a small AT. A more interesting result is that, 
for the negative T-process with T\ = 0.7, T2 = 0.4 signifi- 
cant violation of Eq.(|l.l|) has been observed, while in the 
corresponding positive T-shift protocol Eq. ( [iTl| ) is satis- 
fied within accuracy of the present numerical analysis. 
The former non-cumulative memory effect has been, for 
the first time to our knowledge, observed in simulations 
on the 3D EA models. Although the phenomena appear 
asymmetrically with respect to the direction of temper- 
ature changes, the deviation from the cumulative mem- 
ory scenario in the negative T-shift appears qualitatively 
similarly to the one observed recently in experiments on 
the AgMn spin glassriP We tentatively regard them as 
an effect associated with thcJ^pmperature-chaos effect in 
the equilibrium SG phase, E2rE2P and call thcmJJnaJchaos 
effect'. In contrary, as in the previous workJHJ'EH such 
rejuvenation effects observed experimentally just after 
the T-shift have not been detected even with AT = 0.3 
in the present work. 

The organization of the paper is as follows. In the next 
section we describe our strategy of the simulation and the 
method to evaluate x"(w;i) from the spin autocorrela- 
tion function making use of the fluctuation-dissipation 
theorem. In §3 we explain how to specify the effective 



waiting time from the obtained data of x"( w ; £)> an d 
present the results of t^, or i?T 2 (^D for T-shift pro- 
cesses with various sets of T\ and T2 . In the final section 
we discuss our results, emphasizing on the length and 
time scales involved as well as on the relation to the ex- 
perimental results. 

§2. Method of Analysis 

Numerical simulation on a well-defined microscopic SG 
model, such as the EA model investigated in the present 
work, is of quite importance in studying aging phenom- 
ena. It enables us to faithfully realize any T-shift pro- 
cess and observe any quantity in principle. For example, 
the SG coherence length £,T[t]{t), which we regard as the 
mean SG domain size i?T[t](4L has been calculated from 
the replica overlap function.^ For the isothermal pro- 
cess of the 3D Ising Gaussian EA model which we study 
in thfi—ppBent work, it is well described by the power 
lawEJ'ES'E^ written as 



R T (t)/lo = b T (t/t Q ) 1/z{T \ 



(2.1) 



where Iq and to are microscopic length and time scales 
(lo = 1 lattice distance and to — 1 MC step for simulated 
results), bx is a weakly T-dependent constant, and the 
exponent l/z(T_) linearly depends on T except for the 
region near T C .E5 

The growth law of Rr(t) in isothermal aging different 
from Eq.(2T) was_proposed in the droplet theory due to 
Fisher and HuseO It is written asO 



R T (t)/L Q (T) = L(t/T (T)), 



with the scaling function L(x) given by 

L( X )r 



x 1/z (x «: 1), 

log 1 /^) (x»l). 



(2.2) 



(2.3) 



Here T (T) (~ lot~ u ) is the crossover length, r (T) (~ 
toe~ zu ) is the attempt time for thermal activation pro- 
cess of droplets, where e = (T c — T)/T c , and z and v are 
the critical exponents associated with the criticality at 
T c . The exponent i[> m Eq.(2.3) is, on the other hand, 
intrinsic to the activation dynamics of droplets. We have 
recently confirmed the above growth law, including the 
crossover from the critical dynamics (x <C 1) to the ac- 
tivated dynamics (x < C-1), in the EA SG model but in 4 
dimensioned (see alsoEJ). However, the simulated data 
of Rrit) in the 3D EA model are compatible with both 
the power, law of Eq.(2.1) and the logarithmic law in 
Eq.(2.3)0§> 

Our strategy in the present work is as follows. Be- 
cause of the circumstances of the 3D EA model men- 
tionedjust above as well as those of the recent experi- 
ment ,£3 we do not go into the question which growth law 
is a correct one for Rr(t) observed in the time- window 
( < 10 5 MCS) of t he simulations, and simply use our 
results, i.e., Eq.(2.1), in relating a time scale of observa- 
tion to a length scale of the SG short-range order. Then 
we examine whether the cumulative memory scenario is 
sufficient or not to interpret the obtained results of the 
length scales. We also restrict ourselves to the tempera- 
ture range of T/J = [0.4, 0.7] as already noted in §1. In 
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this temperature range l/z(T) in Eq.( |2.1| ) is well propor- 
tional to T and the dynamics is considered to be domi- 
nated by the activated process. But the prefactor br still 
exhibits weak dependence on T even in this range which 
will turn out not to be neglected in our present analysis. 
Lastly we have examined T-shirt processes with various 
values of the temperature differences AT, in particular, 
a larger one than that studied in our previous work.ES' 

In the aging study at T through the ac susceptibility 
at frequency ui we need to introduce another length scale, 
Lritu,), with = 2tt /ui. It is a mean size of spin clusters 
(or droplets) which can respond to the ac field at T. In 
the droplet picture the aging (or ^dependent) part of 
x"(a->; t) in an isothermal aging is described by a function 
of L T (^)/i? T (i).B0) As will be discussed in §4 below, 
x"{u;t) in some T-shift process to T at t — t w is given 
in terms of Tt(^) / Lt{t) where r = t — t w . Thus it 
provides us information of Rr(t) or Lt(t) in the aging 
process since Lx(t u ) is independent of t or r. 

By the experimental condition of measuring x"(w;£), 
t or t is necessarily larger than t u in general. This 
time regime is called the quasi-equilibrium one, where 
the fluctuation-dissipation theorem 1IT1T) is expected 
to hold well, though approximately.ESEZP Therefore, in 
the present work, the out-of phase component of ac sus- 
ceptibility, is evaluated from the spin auto- 
correlation function 



C{r;t) = (Si(T + t)Si(t)) 
via the FDT aS 

7T 8 



2T dlnr 



C(r;t) 



(2.4) 
(2.5) 



In Eq.(2.4), Si(t) is the sign of the Ising spin at site i 
at time t which is measured in unit of one MC step. 
The over-line denotes the averages over sites and over 
different realizations of interactions (samples), and the 
bracket the average over thermal noises (or different MC 
runs). In this evaluation of x"( u, jt) we are completely 
free from any nonlinear effect of the ac-field amplitude. EP 
In our previous workE^r we studied the susceptibility 
defined by 



X(u;t)=-[1-C(r;t)] 



(2.6) 



It is just the ZFC susceptibility: the induced magneti- 
zation (divided by h) at an elapsed time of r under the 
field h which is switched on after the system has aged 
under h = by a period of t. For slow processes of our 
present interest, x(w;i) ' xs essentially regarded as the in- 
phase component of the ac susceptibility, \'(uj;t). As 
in the experiments, simulated x"(w;i) exhibits larger ef- 
fects of aging relatively to its own absolute magnitude 
than x'(u>;t) or x(w;t) docs. However we have to nu- 



merically evaluate the logarithmic derivative in Eq.(2.5) 
to estimate x"{ UJ 'ii)- bi the present work we have cal- 



culated several C(r;t) in Eq.(2.4), each of which is the 
average over one MC run for each sample but typically 
over 1600 samples. The linear system size is fixed to 
L = 24. The error bars on x"(w;i) drawn in the figures 



shown below indicate the variance in the results of the 
numerical derivative on these several sets of C(r;i). 

§3. Results 

3.1 Isothermal aging 

Before going into discussions on the T-shift protocol, 
let us here present the results of x"(w; t) obtained in the 
isothermal aging. In Fig. 1 we show x"( w ; t) with t u = 64 
in the isothermal aging at various temperatures. They 
play an important role in the following arguments on the 
T-shift protocol, and we call them the reference curve at 
each temperature. 
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Fig. 1. x"(u^: () with tu = 64 in the isothermal aging at various 
temperatures. 



In Fig 2 x"(o>;*) at T = 0.6 for t u = 16 ~ 256 
are shown. The filled symbols are raw data plotted di- 
rectly against t, while the open symbols are the same 
x"{u\t) plotted against u>t with no vertical shifts of the 
data sets. All the sets of data thus plotted nicely col- 
lapse to a universal curve. This reconfirms that the 
wt-scaling—Df x"(a;;t) holds also in the present model 
spin glassEj as observed experimentally .EP As pointed 
out in §2, the time evolution of is considered 

to be a function of Lt(C)/-Rt(*) m the droplet picture. 
The wi-scaling th en comes out from the first equality 
of Eqs.(lT) and (^.l|) S3 The response in equilibrium, 
x"q(w) = Unit— voo x"(^; t), is hardly extracted from our 
present data. We could not detect even its relative dif- 
ference with cu, which should be reflected as the vertical 
shifts of the data sets in the above scaling analysis. 



3.2 T-shift protocol 

In Fig. 3 we show x"(o;;i) with t u = 64 numerically 
observed in the negative (positive) T-shift protocol. The 
temperature is changed from T\ = 0.7 (T2 = 0.5) to 
T2 = 0.5 (Ti = 0.7) at different waiting times t w i (t W 2). 
The observation starts from t = t wi + t w after each T- 
shift. Similarly to t) previously investigated,E£P we 
see the following characteristic features, 
a) Each x"(a;;t) rapidly undershoots (overshoots) the 
T2 (Ti)- reference curve. 
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Fig. 2. x"(tj;t) in the isothermal aging at T = 0.6. The filled 
symbols are the data plotted directly against t, while the open 
symbols are the same data plotted against dAt/t^. 



b) merges to the T 2 (Ti preference curve from 

below (above). 
For the negative T-shift, in particular, the value of 
x"{w,t) just after the negative T-shift is the larger rela- 
tively to the T 2 -reference curve, the larger is t w %. How- 
ever, ~x"(Lo\t) just after the shift does not exhibit over- 
shooting of the Ti-reference curve, a phenomenon which 
we call the strong rejuvenation effect in the present pa- 
per. This is also the case even if we .dare to measure 
t) at r = t — t w i smaller than t u M3 



T-shift protocols 
r 1 =0.7,r2=0.5,-£O-'=64 



r wl =512 >~e- 

r„i=2048 > *-■ 
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value of tshi (— 3300 — 512 in the figure) it merges to 
the reference curve and lies on it afterwards. We regard 
time t required for the shifted branch to merge to the 
reference curve in this situation as the effective waiting 
time, t^, introduced in §1. If £ s hi is further increased 
the branch merges to the reference curve from above and 
again at larger t than t^. Thus the chosen t s hi which 
corresponds to t^ yields the shortest time for the shifted 
branch to merge to the reference curve. 

An interesting observation in the above analysis is 
that the time at which the properly shifted branch 
merges to the T2-reference curve is nearly equal to 2f^; 

*wi + W + Ci — 2i wi and so *wi - *wi + *shi- This 
aspect, which is by no means trivial, has been commonly 
observed in most of the T-shifts examined in the present 
work. Using this observation, we estimate errors of t™[ 
as follows. We judge by eyes the largest i s hi for which 
the shifted branch of x"(w; t) certainly crosses with but 
not merges to the reference curve, and this value of t B hi 
gives a smallest estimate of tf^. Similarly the smallest 
i s hi for which the shifted branch merges to the refer- 
ence curve at t > 3tf^[ yields a largest estimate of . 
The three shifted branches shown in Fig. 4 correspond 
to these smallest, mean, and largest estimates for t^[. 
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Fig. 4. x"(u);t) in the negative T-shift protocol from T\ = 0.7 
to T2 = 0.5. Each symbol represents the branch of \"(uj\i) at 
t > t w i shifted by an amount of t s ^i indicated in the figure, the 
line with the smaller symbols is the T2-reference curve. 



Fig. 3. x"( UJ it) in aging with negative and positive T-shifts be- 
tween Ti =0.7 and T2 = 0.5 at t = indicated in the figure. 
The upper and lower curves with the smaller symbols represent 
the Ti- and T2-reference (isothermal) curves, respectively. 



3.3 Effective waiting time 

Feature b) above is examined in more details in Figs. 4 
and 5. We note that the i-axis in these figures is linear 
in t. Within the time window of Fig. 4 the merging of 
bare x"i UJ it) (denoted by t s hi = 0) to the T 2 -reference 
curve is not seen. If, however, the branch of x"(w;t) 
at t > t w i is shifted to the right by an amount denoted 
by i s hi, it crosses the reference curve and merges to it 
at a smaller t than that with t s hi =0. At a certain 



The above analysis for the negative T-shift protocol 
also works for the positive T-shift protocol. A typical 
example from T 2 = 0.5 to T\ = 0.7 with t w2 = 16384 is 
shown in Fig. 5, for which we have to reread the suffix 1 
by 2 and vice versa in the above argument. Also in this 
case the branch of x"(a>;t) at t > i w2 is shifted to the 
left by tgitf. A too large £ s h 2 makes the shifted branch to 
overshoot the Ti-reference curve, while a too small i s h 2 
significantly delays the merging. With a proper chosen 
t s h2 we obtain t^ 2 (— 1600 from the figure) which satisfies 
t^ 2 = tw2 ~ tsh2 ■ Its error bar is similarly evaluated from 
the two other t s h 2 's indicated in the figure. 
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Fig. 5. x"( w i t) with tu> = 64 in the positive T-shift protocol from 
T 2 = 0.5 to Ti = 0.7 with t w2 = 16384. The data in a large time 
scale are shown in Fig. 3. The three sets of symbols represent 
the branches of x"( ur >t) at i > i w 2 shifted to the left by the 
amount t a h2 indicated in the figure. The line is the Ti-reference 
curve. 



AT = 0.1 the effect of the T-dependence of 6t is neg- 
ligibly small. The effect is, however, significant for 
AT > 0.2. Thus the T-dependence of &t has to be prop- 
erly taken into account to judge the cumulative nature 
of memory observed even in the temperature range ex- 
amined in the present simulation. 
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3.4 Cumulative memory and 'chaos' effects 

In Fig. 6 we plot (t W 2) as a function of t w i (t^|) 
obtained in the previous subsection in the negative (pos- 
itive) T-shift protocol for three sets of (T 1; T 2 ). Here we 
have followed the idea of 'twin-experiments' in the re- 
cent work.EH' Before the explanation of the lines drawn 
in the figure, we note that the data points of both neg- 
ative and positive T-shifts with T\ = 0.7, T 2 = 0.6 are 
seen to lie on a certain common curve, while this is not 
the case for those with T x = 0.7 and T 2 = 0.4. The for- 
mer is expected from the cumulative memory scenario. 
But the latter data points clearly indicate a violation to 
the scenario irrespectively of the growth law of the SG 
domains. 

Now let us explain the lines in Fig. 6. The solid ones 
represent the relation between the (t W 2) and i w i (t^|) 
when the cumulative memory scenario represented by 
Eq.(Ll) (the one whose suffix 1, 2 interchanged) com- 
bined with the growth law of Eq.( |2.l|) holds. For the lat- 
ter we have explicitly used the following sets of the pa- 
rameter values (T;z(T),6p) we previously obtained:t3 
(0.7; 8.71, 0.779), (0.6; 9.84, 0.782), (0.5; 11.76, 0.800) 
and (0.4; 14.80, 0.818). As seen in the figure, for a small 
AT (= 0.1), both t^ and ff 2 lie on the solid line. The re- 
sults confirm the cumulative memory scenario described 
in §1. With AT = 0.2, tf 2 still lie on the solid line 
but tends to deviate, though a little, from it. For 
AT = 0.3, significantly deviate from the solid line, 
which is incompatible with the cumulative memory sce- 
nario. It should be emphasized, however, that the cor- 
responding x"(w;i) (not shown) in this process exhibit 
features a), b) mentioned in §3.2, i.e., no strong rejuve- 
nation. 

If the wea k T- dependence of &t is discarded the con- 
dition of Eq.(l.l) is reduced to 

Ti/T 2 

(3.1) 



L wl 
TO 



t w 1 
TO 



Fig. 6. vs t w ; for the negative (i = 1) and positive (i = 2) 

T-shift protocols for three sets of Ti and T2 indicated in the 
figur e. T he lines represent the expected behavior from Eqs.(l.l) 
and (2.1) as explained in the text. 
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which is also shown by the dotted lines in Fig. 6. For 



Fig. 7. Relations t*5 vs. t w j drawn by means of t he c orrespond- 
ing domain sizes evaluated by Eq. ( |l.l| ) with Eq. (2.1). Here we 
omit the error bars. For all the sets of Ti,T 2 their magnitudes 
are comparable with those shown in Fig. 6. 



In Fig. 7 we replot our data in Fig. 6 as well as those 
of other sets of (T 1; T 2 ) in terms of the lengths, where 
Ri = Rt ± (i w i ) and Rl s = Rt 2 (t^{ ) which are evaluated 
by Eq.(2.1) using t w i and extracted in the negative 
T-shift processes (i? 2 and i? 2 ff in the positive T-shift pro- 
cess are similarly evaluated). The line in the figure is 
what is expected from the cumulative memory scenario, 
i.e., i?° ff = Ri for both i — 1,2. We see clearly that 
this is the case for both negative and positive T-shift 
processes with AT = 0.1 within the time window of the 
present simulation. For AT = 0.2 the data of the pos- 
itive T-shift satisfy the condition i? 2 ff = i? 2 , but those 
of the negative T-shift exhibit the tendency i?J ff < R-y. 
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Behavior of the T-shift with AT = 0.3, i.e., Zi = 0.7 and 
T 2 = 0.4 is as already described above and is interpreted 
below to be due to the 'chaos effect.' 

According_tn_the theory for the temperature-chaos in 
spin glasses,OEiP the SG equilibrium configurations at 
different temperatures, Xi and T 2 , are completely uncor- 
rected with each other in the length scale larger than 
I at °c AT -1 / 1 ', where I at is called the overlap length 
and £ (> 0) the chaos exponent. An important prob- 
lem here is how the existence of Iat affects the non- 
equilibrium aging dynamics. Let us consider a negative 
T-shift process with AT, for which Iat is supposed to 
be sufficiently smaller than i?Ti(iwi), and introduce the 
time scale t ov i by Tt 2 (^ovi) = Iat- At a time range after 
the T-shift specified as tf^ S> r ^> t ov i, a longer part 
of the memory imprinted before the T-shift is still pre- 
served, but such a memory is expected to be irrelevant 
to the equilibration process to the SG ordered state at 
T = T2. Thus the system looks as if it is already in 
the isothermal aging state at T 2 . Then, if our analysis 
to determine t^i is applied to this T-shift process, the 
expected result is ~ t ov l irrespectively of t w i . The 
circumstances are the same for the positive T-shift proto- 
col. Consequently, R° in Fig. 7 is expected to saturate 
to Iat at large Ri both for i = 1,2, and the data for the 
negative and positive T-shifts come out symmetrically 
with respect to the line of i?® ff = Ri . 

We tentatively interpret our data of the negative T- 
shift with AT = 0.3 as an early stage of the saturation 
described above. Unfortunately, the data are so limited 
that we cannot figure out a value of Iat- Also the corre- 
sponding positive T-shift data nearly coincide with the 
line R° s = Ri, i.e., the two sets of data are by no means 
symmetric with respect to the line. One of the reason of 
this asymmetric behavior may be due to our method to 
specify the effective aging time combined with the time 
scales in our simulation. Although Iat is common to 
the negative and positive T-shifts, the separation of time 
scales i w2 , i v2 and £^ 2 i n the positive T-shift is much 
smaller than that of t w i , i ov i and in the negative T- 
shift. This is due to a large difference in the growth rates 
at the two temperatures. It is then rather hard to detect 
a possible small difference between i ov2 and within 
our present analysis. With these reservations, we inter- 
pret our results of the T-shift process with AT = 0.3 as 
a dynamic process which reflects the temperature-chaos 
predicted for the equilibrium SG phase. 

§4. Discussions 

In the T-shift protocol examined in the present study, 
the cumulative memory scenario has been confirmed for 
T-shift processes with a small magnitude of the shift, 
i.e., AT = T x - T 2 = 0.1. This has been done by close 
comparisons of \"{u>;t) after the T-shift with that in 
the isothermal aging at T2 (reference curve). However, 
there have been little experiments^hich directly mea- 
sure tf^y similarly to our_analysisc3 : E£P An example is 
the one by Mamiya et alf^l who analyzed the aging dy- 
namics in the SG-like phase of a ferromagnetic fine par- 
ticles system. In the T-shift process with Ti = 49K and 
T 2 = 47K (with T g ~ 70K) they observed ^ - 3 x t wl 



for t w i = 2.0, 15.0ks. If we suppose t = 10 6 $-£or a 
magnetic moment carried on by each fine particle,EU we 
obtain ~ (2.4 - 2.7) x < wl from Eq. ([[_!]). The result 
is rather satisfactory and implies that the cumulative 
memory scenario works as well for the T-shift process 
with a small AT in this SG material. 

Our results on the negative T-shift protocol with 
AT = 0.3 have turned out to be incompatible with the 
cumulative memory scenario. The period necessary 
for the system to become in a T 2 -isothermal aging state 
after the T-shift is signific ant ly smaller than the valu e 
of estimated from Eq.(Ll) combined with Eq.(2.1). 
In this negative T-shift process which violates the cu- 
mulative memory scenario, however, the strong rejuve- 
nation phenomenon just after the T-shift, which is de- 
scribed in §3.2, has not been detected. We have therefore 
attributed the non-cumulative memory effect we have 
found to the 'chaos effect'. 

Quite recently Jonsson et al. (JYN) have reported 
the chaos effect which symmetrically appears in the pos- 
itive and, negative T-shifts in a Heisenberg-like spin glass 
AgMnJiiJ They have measured the ZFC magnetization 
with schedules of temperature change corresponding to 
the T-shift protocol discussed in the present work but 
within a very small range of AT (< 0.012T C ). The 
logarithmic-time derivative of the ZFC magnetization, 
S(t;t w i), exhibits a peak, whose position is considered 



to be at r 



t - t 



wl 



cir 



the time required for the 



merging to an isothermal state at the new temperature 
just investigated in the present work (see the discussion 
below). In fact, our Rf -vs-Ri plot in Fig. 7 of the neg- 
ative T-shift with AT = 0.3 is in qualitative agreement 
with their LgS-vs-Lxiitv,) plot, where their L e g (Tti (^w)) 
just corresponds to our i?J ff (Ri)- In contrast to our nu- 
merical observation, however, their data for the positive 
T-shift plotted in our way appear symmetrically to the 
negative one with respect to the line i?f s = Ri . The over- 
lap length Iat estimated by scaling analysis has turned 
out to be larger than Ri (= Lt^w)), or before the sat- 
uration of R° mentioned in §3.4. One of the reasons of 
this discrepancy between their experiment and our simu- 
lation may be the Heisenberg-spin nature in their mate- 
rial AgMn; Heisenbjecg spin glasses are more chaotic than 
Ising spin glasses.li§E!P This point is very interesting and 
to be further pursued. 

Next let us make a few comments on aging at a time 
range r = t — t w i < after the T-shift, which we_have 
called the transient regime of the T-shift process.liS A 
main idea for this regime is the droplets-in-domain sce- 
nario which involves at least two characteristic length 
scales as mentioned in § 1. One is the mean domain 
size at r — 0, i.e., i?Ti(iwi) and the other is Lt 2 (t), 
the mean size of droplets (or subdomains) which are al- 
ready in local equilibrium of the shifted temperature T 2 
at time r after the T-shift. Associated with the growth 
of Lt 2 (t) some peculiar features have been observed. An 
example is a non-monotonic time evolution of the energy 
density in a positive T-shift process (see Fig. 4 inH). 



It is recently named as the Kovacs effect inE!P since the 
qualitatively similar phenomenon was first observed in 
polymer glassesO 
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In the droplets-in-domain scenario here we implicitly 
assume that droplets (or subdomains) in local equilib- 
rium of T2 distinguish themselves from those in local 
equilibrium of T\, On the other hand, the Kovacs ef- 
fect has been observed in negative T-shift processes with 
AT = 0.2, for which the 'chaos effect' is not clearly 
detected in Fig. 7, or 1& T > i? Tl (i wl ). This strongly 
suggests that in nonequilibrium aging dynamics the spin 
configurations which we have so far supposed to be in lo- 
cal equilibrium at two different temperatures differ from 
each other even in length scales shorter than the equilib- 
rium I at of the corresponding temperatures. This view- 
point is in contrast to the argument of the temperature- 
chaos in equilibrium, and is worthy to be investigated. 

From the droplets-in-domain scenario mentioned just 
above, the strong rejuvenation experimentally observed 
in the ac susceptibility measurement discussed in §3.2 
can be regarded as one of such peculiar phenomena in 
the transient regime. As mentioned in §2, the ac re- 
sponse associates another short length scale Lt 2 I n 
the time range t u <C r <C t^j so that Lt 2 (t u ) < Lt 2 (t) < 
-RTi(iwi) holds, the t-dependent part of x"(w;t) is gov- 
erned dominantly by the ratio LT 2 (t u ) / Lt 2 (t) since 
droplets responding to the ac field less feel the existence 
of larger domains of i2r 1 (f w i) than that of smaller sub- 
domains of Lt 2 (t). The strong rejuvenation alone is 
therefore not necessarily incompatible with the r.iimu.- 
lative memory scenario. In neither the previoual3'l23' EiP 
nor the present simulations, however, the strong rejuve- 
nation has been detected. This may be again attributed 
to the smallness of the time window; the separation of 



t u < T <C tl 



is not enough in the 



the time scales, 
simulations. 

Quite recently Yoshino and the present authors have 
argued based on the numerical results on the 4D Ising EA 
model that fluctuations of droplets, whose size becomes 
comparable to that of the preexisting domains, become 
anomalously large, and that they are responsible to the 
occurrence of a peak in S(t;t w ) of the isothermal ZFC 
magnetization at r ~ t w where r = t — t w is_the time 
elapsed after the measuring field is applied.c.1 At the 
end of the transient regime of the T-shift process, or at 
the merging to a T2-isothermal aging state, similar large 
fluctuations and so a peak in S(t;twi) at r ~ tf~[ are 
expected to appear so long as AT is relatfpJy small. In 
fact this has been experimentally observedt3 as already 
mentioned above. 

Combining the arguments based on our numerical re- 
sults, in particular, the droplets-in-domain scenario and 
the cumulative memory one, with a possible existence of 
the 'chaos effect', we can think of the following behavior 
that a spin glass exhibits in the negative T-shift protocol 
of aging depending on the magnitude of AT (and simi- 
lar behavior also for the positive T-shift protocol) . For a 
sufficiently small AT, the merging to the T2-isothermal 
aging state is observed at r ~ V^ x , where is given by 



Eq.(Ll), i.e., by the cumulative memory scenario. When 
AT becomes large, both the strong rejuvenation just af- 
ter the T-shift and the merging to the Tij-isothermal ag- 
ing state at r ~ are expected to be observed. The lat- 
ter, however, becomes to be hardly detected by such an 



ac susceptibility analysis done in the present work since 
LT 2 (t w ) is much smaller than Lx 2 (^S)- Instead, it is ob- 
served, through a peak in S(t; t w i) of the ZFC magnetiza- 
tionEif The extracted value of in this case either sat- 
isfies the cumulative memory scenario (iovi ^wi ) or ^ 
already strongly affected by the chaos effect (t ov i ~ ) . 

What is the expected behavior for T-shift processes 
with a sufficiently large AT for which i ov i <C or even 
iovi < C b in holds? Here t^ n is the shortest time that 
the temperature T2 is experimentally stabilized after the 
T-shift. The recent experimental results are claimed to 
reach this regimeJlj'H and a theory for the chaos effect 
on T-shift and T-cycling processes in this regime has 
been proposed by Yoshino et alS3 Unfortunately this 
regime has not been realized in the numerical simulations 
on the 3D EA model. Probably it needs a sufficiently 
large t w \ , larger than the time- window of our simulation 
( < 10 5 to), to realize the condition i?Ti(*wi) > 'at- In 
order to further explore aging dynamics in the T-shift 
protocol, one has to systematically choose values of the 
parameters AT and t W i even in experiments, since their 
time-window is similarly small (~ 5 decades) to that of 
the numerical simulation though its absolute magnitude 
is large (Is ~ 10 12 i ). 

To conclude, we have numerically studied the T-shift 
protocol of aging in the 3D EA spin-glass model through 
the measurement of the ac susceptibility. For processes 
with a small magnitude of the T-shift, AT, the memory 
imprinted in the first stage of isothermal aging is pre- 
served under the T-shift and the SG short-range order 
continuously grows with a rate intrinsic to the tempera- 
ture changed (cumulative memory scenario) . For T-shift 
processes with a large AT the deviation from the cumu- 
lative memory scenario has been observed for the first 
time in the numerical simulation. We attribute the phe- 
nomenon to the 'chaos effect' which, we argue, is qualita- 
tively different from the so-called the rejuvenation effect 
observed just after the T-shift. 
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